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ABSTRACT 


A modern area of fluid dynamics, called nanofluid dynamics has 
become increasingly popular in novel hematological treatments. 
Motivated by these recent developments in the field of 
nanotechnology, a theoretical model for catheterized stenosed artery 
having permeable walls, treating blood as a base fluid for various 
shapes of nanoparticles suspended, is presented. The disposition of 
blood is described to be that of a viscous nanofluid. The calculations 
are executed using mild stenosis conditions with heat and mass 
transport phenomenon. The varying shapes of nanoparticles which 
include bricks, platelets, cylinders and blades have been given due 
attention using Hamilton-Crosser equation that used to obtain 
expressions for thermodynamic properties like thermal conductivity 
and viscosity using the appropriate boundary conditions. The 
comparative effects for velocity and temperature have been depicted 
graphically. This physical model has useful application in nano-drug 
delivery for the cure of various cardiovascular diseases. 
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1. INTRODUCTION 

The blood motion analysis through unhealthy arteries 
is a crucial domain of bioengineering research and 
analysis. Enhanced clinical data for mathematical 
models and recent developments in the computational 
tools has attracted ample consideration in the 
contemporary years. The most recurring arterial 
ailment is atherosclerosis. The unusual or unhealthy 
development in wall of artery at distinct locations in 
circulatory system, in medical terms is defined as 
stenosis[1]. Apparently, the enlargement of stenosis 
in an artery lessens in flow of blood. The growth of 
such plaques in blood vessels can cause sporadic 
blood flow rate, high arterial wall shear stress and 
boundary layer detachment[2]. Ku [3] in his review 
mentioned about the energy loses, turbulence and 
other aspects of hemodynamics. 


Numerous assessments established on the principle of 
dynamics of fluid address the problems on blood 
motion in composite stenosed artery having a 
catheter. Catheter is designed using medical grade 
polyvinyl chloride and _ polyester constructed 
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thermoplastic polyurethane etc. [6].The incorporation 
of catheter effects the hemodynamics. Srivastav [9] 
explored composite stenosis in presence of a catheter 
with permeable walls. Srivastava and Srivastava [10] 
analyzed a stenosed artery with catheter as well as 
without catheter. Srivastava et al [11] studied the 
problem of composite stenosed artery with 
catheterization considering blood modeled as non- 
Newtonian fluid. Ellahi et al [12] explained flow of 
bloodin arteries with composite stenosis while 
characterizing blood as a micropolar fluid. 


The fast pace technological evolution over the past 
years has caused increasing stipulation for novel 
medical treatment methods, surpassing — the 
conventional ones, to improve the therapeutic effects 
and lessen the side effects. Nobel laureate Richard P. 
Feynman introduced the term nanotechnology, and 
since then there have been various revolutions in this 
field[13]. Miscellaneous applications in energetics 
and bio-medical science of nanoparticles in a base 
fluid has been classified under nanofluid dynamics, a 
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new branch of fluid dynamics. The term metal 
nanoparticles describe nano-sized metals with 
dimensions within 1-100 nm[14].Nanoparticles are 
used as new therapies for the treatment of 
cardiovascular diseases due to high target 
specificity[15]. Choi [16] explored the nanofluids. 
Ellahi et al [21] studied nanoparticle outcomes in 
arteries having composite stenosis with permeable 
walls. Nadeem and Ijaz [22] theoretically examined 
nanoparticles carrying drug in stenosed artery 
accompanied by slip effects on the wall. Chatterjee et 
al [23] used Bernstein polynomial estimation to study 
the effects of nanofluids in stenosed artery. Rathore 
and Srikanth [24] considered blood as a micro-polar 
nanofluid in an artery including stenosis under in the 
influence of catheter. 


Since the chief channel of heat transmission in fluids 
is convection, its efficiency relies on the thermo- 
physical attributes of the fluids[25]. Nanoparticles 
possess very high thermal conductivities. There are 
diverse geometries associated with them. Various 


2. Formulation 


particle geometries either occur commonly or are 
schemed for distinct uses. Particle configuration is a 
significant nanoparticle criterion. They critically 
effect the nanofluid performance. In much research, 
spherical-shaped nanoparticles have been considered. 
However, practically they have limited applications. 
Drug delivery or clinical diagnosis have a usage of 
non-spherical shaped nanoparticles. Evidently, 
nanoparticles of geometries like brick, platelet, blade 
and cylinderetc have colossal relevance in drug 
delivery[27]. Timofeeva et al [28] studied thermal 
conductivity and viscosity of varying geometries of 
nanopatticles. 


The current article is devoted to analyze the 
consequence of varying shapes of nanoparticlesin an 
artery having composite stenosis with a catheter. This 
is a novel comprehensive examination of mass and 
heat transfer of nanofluid flow in an artery having 
composite stenosis artery in presence of catheter with 
due attention given to various shapes of nanoparticles 
and various thermodynamic properties of nanofluids. 


Steady and laminar flow of blood in a cylindrical artery of length L has been considered. The incompressible 
flow in the artery is modeled using cylindrical co-ordinates. Let the velocity vector be (u', v',w") here u' is inz‘ 
direction, 1" is along radial direction orr‘direction andw'is along é" direction. For the case of axis-symmetricity 
w' = Q, since 6" represents the circumferential direction. The configuration of composite stenosis is defined as: - 


Ro —! macs '—d) rdsz'sd aes 
R'(z") =4Ro- ze +cos— = (z" —ad—)) ;d+Bsz'sdt Ro (1) 
Bo ; otherwise 


where @ represents location of stenosis, £, is length of the stenosed part, Lis arterial length, Ry is radius of 
normal artery and R'(z‘) is arterial radius with stenosis.R , is catheter’s radius (Figure 1) 


Figure 1 
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The governing equations with regard to nanoparticles are given as: - 


454 =0 (2) 
i ("5+ eo) ~ a lng Gore +33 a a (3) 
Put ('& ill e) = o oe +t Png ( * o = =.) + 9Pnf¥ng(T' — To) oe 
@w'S+wS= —t_ Gat $2 ++ (5) 


Pf Bp, Paf Png 


T,, is temperature of wall of catheter and T, is temperature of artery. T' is nanofluid temperature. The heat 
generation or heat absorption parameter is represented by H. Now considering nanofluid model, ,,- represents 
nanofluid viscosity, k,,,- is nanofluid thermal conductivity, ,,- is nanofluid density, g,,-y,,¢ 1s nanofluid thermal 
expansion coefficient ae Pal, - is nanofluid heat capacitance. 


Boundary conditions are specified as: - 


uw =u',,T' =T,at rr’ = R(z'} (6) 
u =0, T'=T, atr' =R, (7) 
aur og a ' ,_ wire 

art = Jeg a — Up Jat = R"(z") (8) 


u', is velocity of blood at the wall of catheter and =. is slip velocity of arterial boundary that is calculated using 
Darcy law. is a dimensionless parameter depending upon the nanofluid and artery. Equation (8) is obtained 
from the Beavers and Joseph condition [35] for boundary between a free fluid and a porous medium. 


== ky, dp! <j a 
uy —= re, + GP ag Yang qr" r Ty) (9) 


Using T' =T gat r' = R"(z") 
a. (10) 


Despite extensive research, there are not many reliant universal theoretical models to define the nanoparticle 
thermal conductivity. Maxwell introduced the macroscopic Effective Medium Theory(EMT) which was further 
developed by Hamilton and Crosser for aspheric particles. This model states the thermal conductivity of two 
components in a mixture which is an outcome of conductivity of unblended materials, their respective 
compositions and the fashion in which they are spread in the mixture. It calculates the thermal conductivity 
(k,, ) of nanofluid. It comprises empirical shape factor m = 3,/1, where fis the sphericity which is termed as 
the ratio of the sphere’s surface area to the real particle’s surface area having equal volumes. The values of n for 
various shapes of nanoparticles is listed in table 2. 


Kng _ tla-Dikgtin—-Liky—ke eh (11) 
ke ~ kept (n—Dke—-(Rp—-k pe 


where kis the thermal conductivity of blood, while k,,is nanoparticlethermal conductivity and¢ is nanoparticle 
volume fraction. 
The nanofluid density is expressed as 

7 =A-4)e, + dp, (12) 
where Pr is density of blood and Prey is nanoparticle density.Similarly, nanofluid specific heat capacity is written 
as 
Parag (1—@)p,e, e + PP py (13) 


The nanofluid thermal expansion is measured as 
Png ¥ng = (1— P)epye + PPy¥ (14) 
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where Prop, and Pr¥pis the specific heat capacity and thermal expansion of blood. Pye, and Pa¥pis 
nanoparticle specific heat capacity and nanoparticle thermal expansion. 

The nanofluid viscosity is given as 

Eng = (14+ A@) py (15) 

where Ais specific to the shape of nanoparticles in the fluid(Table 2). ,,- is nanofluid viscosity and pr, is blood 
viscosity. 


The non-dimensional variables are:- 


r F r e fq Rotgpg # avg egR o(T,-Ty) 
r=—,z=—,v= ¥ u= u p= RgP! o=-, R, Rotavg OF 6.= fer ote =. 
Ry Ry Varg Varg Fang Ry MF ip HF 
a HR r'—t, Ie 
é=1,,=—s_ g=7-3, pa=*2 
0 (T-Th) ke Tet Ay (16) 


In equations (16), &, is Grashof number, uw, is average velocity, R, is Reynolds number ,f is heat source 
parameter andDa@isDarcy number. The blood current in arteries with small diameters is sluggish therefore 
intensity of inertial forces is negligible hence there remains only one component parallel to the axis, as a result 
v' = 0. Using the mild stenosis condition d = 6" /R, << 1, the transformed equations (2) to (5)are stated as:- 


z= wm 
geet = FS +iF +a $) + o2RI(1+ A8)6,0 (18) 


ae 186 keytln—L) kp-(ky—k ‘f (19) 


kept (n—1)k p+ (a1) (Kp Kk le 


The geometry of the stenosis and the boundary conditions under non-dimensional scheme are:- 


1 —26(z-a); gSzSa+5 
R(z)= 1-2 (1+cos2m(z — @—>)); g+=sz sol (20) 
1; otherwise 
u=u,,8@= 0 atr= R(z) (21) 
u=0,0=1atr=R,/R, (22) 
== = z —(u,—u,) at r = R(z) (23) 
3. Solution 


The solution to the equations (17) to (19) using (20) to (23) is obtained analytically as: - 
The temperature is given as:- 


a 1 kyptln—Dke—(ky—ke lb ja _ 
~ Ginma)-In g/m) gt (n—i)key Hn—D) (yee, pF (z)"(nr 


Ink, /Ry) +r? (-InR(z) + Ink, /Ry) + (n&(z) — Inr)(R,/%y)*)] (24) 


The velocity is given as:- 


(4In R(z) —InR,/Ry) +h 
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_ op (1+Ad) 
dz 4(InFA(s)—ln Fp! Ry) . 
R.2(—-InR(z) + nr)) + (1 + 4g) [a —¢)+¢ sete] G.{ 16(R(z)* (in r- n=) + 
r°(1+In R(z)—-In r)(—In R(z)4+In RB, /R,)+ (In R(z)—In r)(14+In Rez) — 
In R,/Rp) =) [Ses (r*(In R(z))? + 
(R(z)* —r*) In Rk, /R,)(4R(z)* + (3R(2)* — r*) In BR, /R,) + R(z)}? In r( —4R(z)* + 
(-3R(z)*+ 4r7) In 2) + nR(z)(4R(z)*r? + 
R(z)*(3R(z)* — 4r7) In . + (-3R(z)*+ 4R(z)*r* — 2r*) In R,/R,) — 4(r7(InR(z))* - 
2R(z)*In r + (R(z)*- or? +r? Inr )nR,/R, +n R(z)(R(z)* 4+ 1r* — 7? (Inr + 
InR,/Ro)))(R./Ro) +  (InR(z) —Inr) (4+ 3InR(z) —3InR,/Ry)(R,/Ro)*)) — 


(-lnrtlnkg shy hug 
(In&(s)—Ink,/R,) (25) 


(r*(InR(z) —InR,/R,) + R(z)*(—Inr + nR,/B,)+ 


The volumetric flow rate Q is defined as:- 
g= E 2) seru dr (26) 


= dp q,(=) 
Q am dz (legR(s)—legR,)*-¢,(=) (27) 
4. Graphical results and discussions 
Our novel analysis for the flow of various shapes of nanoparticles in a catheterized composite stenosed artery 
considers various parameters like shape factor n,volume fraction ¢, heat source parameter h, Grashof number 
Gr, thermal conductivity k,,, and viscosity y,, ». The results are obtained analytically, and graphs are obtained 
for temperature and velocity. Importantly influence of various shapes of nanoparticles are obtained for cylinders, 
platelets, bricks and blades on temperature and velocity. Fig 2-5depicts temperature of nanofluid versus radius. 
Fig 6-10exhibitsvelocity of nanofluidversus radius. 


Figure 2 shows the temperature @ of the nanofluid in the stenosed artery versus radial direction r for distinct 
values of shape parameters n. The value of shape parameter defines the ratio of real particle’s surface are to 
sphere’s surface area with equal volumes. Greater the shape factor greater the interaction of the given shape of 
nanoparticle with the fluid interface, greater the viscosity. Thus, blade shaped nanoparticles show least viscosity 
because of greatest shape factor. And brick shaped nanoparticles have least shape factor on account of their shear 
thinning behavior with temperature. Consequently, brick shaped nanoparticles show maximum rise in 
temperature while blade shaped nanoparticles show minimum rise in temperature. 


1.2 


—— bricks 
= = = cylinders 
platelets 


seccccece blades 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


Fig 2 Temperature 9 versus radial direction r various shapes of 
nanoparticles 
Figure 3 depicts the variation of temperature 0 of the nanofluid in the stenosed artery versus radial direction r for 
distinct values of volume fractiong for blade shaped nanoparticles in blood. With the enhancement in volume 
fraction for nanoparticles, the number of nanoparticles in blood increases. The increase in number of 
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nanoparticles causes a greater increase in the transportation of nanoparticles via conduction from catheter wall to 
the boundary. Thus, with rise in volume fraction of nanoparticles, temperature of nanofluid rises. 


Fig 3 Temperature 8 versus radial direction r for distinct values 
for volume fraction 


Figure 4 depicts temperature 8 of the nanofluid in stenosed artery versus radial direction r for various values of 
heat source parameter fh for blade shaped nanoparticles in the blood. As the value of heat source parameter 
intensifies it causes an enhancement in heat generation in blood which enhances the value of temperature of 
nanofluid. Thus, temperature profile shows a rise with rise in the value of heat source parameter. 


1.4 
12 

h=4.0 

— O_o ee neces h=5.0 

0.8 h=6.0 

@ 
0.6 
0.4 
0.2 
0 
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 iE 


r 
Fig 4 Temperature 9 versus radial direction r for distinct values of 


heat source parameter h 


Figure 5 shows the graph of temperature & of the nanofluid in the stenosed artery versus thermal conductivity 
k,,-for distinct shapes of nanoparticles. The graph shows that brick shaped nanoparticles have greatest 
conductivity while blade shaped nanoparticles have least thermal conductivity in blood. Thermal conductivity 
enhances with rise in temperature. Thermal conductivity in return depends on shape factor of nanoparticles. The 
larger shape factor causes lesser conductivity. 
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1.2 
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Fig 5 Temperature 9 versuskeffermal conductivity knf for distinct 
values of shape parameters n 


Figure 6 exhibits velocity profile 4 versus radial direction r for various shapes of nanoparticles. Brick shaped 
nanoparticles in blood show highest velocity profile while blade shaped nanoparticles show a least velocity 
profile. The consequence of the shape of nanoparticles on the velocity is because of the viscosity dependence- 
relation of the shape of the respective nanoparticle at a given temperature. Platelets and cylinders have almost 
same viscosity in nanofluid due to elongated structures, thus they show an overlapping profile for velocity. 
Blade shaped nanoparticles show maximum viscosity while brick shaped nanoparticles have least viscosity. 


=—-= blades 


pemeeeeseeeeeettens | 


platelets 
cylinders 


seeseeeee bricks 


1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


r 
Fig 6 Velocity u versus radial direction r for various shapes of 
nanoparticles 


Figure7 depicts the variation of velocity u of the nanofluid versus radial direction rusing distinct values of 
volume fraction @ of blade shaped nanoparticles in the blood. It is seen that the velocity decreases with rise in 
the value of volume fraction of nanoparticles in the nanofluid. This is so because as the volume fraction 
increases, the number of nanoparticles in the blood increases, which makes the nanofluid more viscous. The 
enhancement in viscosity causes an enhancement in the friction force which reduction in velocity. 
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Fig 7 Velocity u versus radial direction r for various values of 
volume fraction } 


Figure8 shows the variation of velocity u of the nanofluid versus radial direction r for distinct values of Grashof 
number (Gr) for blade shaped nanoparticles in the blood. The results show that velocity rises with the rise in the 
value of Grashof number. Grashof number is the ratio of upthrust by fluid density because of temperature 
difference torestraining forces because of viscosity of nanofluid. The increase in value of Grashof number 
reduces the viscosity of the nanofluid which thus increases the velocity. 


ate PS SE oS _ Gr=1.0 
: ee Gr=3.0 
~ Gr=5.0 


-0.5 0/1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 | 


Fig 8 Velocity u versus radial direction r for distinct values of 
Grashof number Gr 


Figure9 exhibits the variation of velocity u of the nanofluid versus radial direction rusing distinguished values 
of heat source parameter fi for blade shaped nanoparticles. The enhancement in heat source parameter causes an 
increase in the temperature which decreases the viscosity of nanofluids hence increasing its velocity. 
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h=6.0 


r 
Fig 9 Velocity u against radial direction r for distinct values of heat 
source parameter h 
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Fig 10 Velocity u versus viscosity nr of nanofluid for various 
shapes of nanoparticles 


Figure 10shows the graph of velocity u against viscosity Bag of nanofluid for various shapes of nanoparticles. 
Brick shaped nanoparticles show highest velocity profile and blade shaped nanoparticles show least velocity 
profile on account of their varying shape factors. The blade shaped nanoparticles show highest viscositywhile 
brick shaped nanoparticles show least viscosity. The increase inviscosity of nanofluid causes a decrease in 
velocity of nanofluid. 


5. Conclusions 

The flow of a nanofluid is modeled in acomposite stenosed artery having a catheter with prime focus on various 
shapes of nanoparticles in blood and various thermodynamic properties of nanofluids. The results obtained can 
be summarized as:- 


1. The nanoparticles having the shape of blades show least rise in temperature while brick shaped nanoparticles 
show maximum rise in temperature. 


The rise is the number of nanoparticles in the blood causes a rise in temperature of nanofluid. 
The increase in heat source parameter causes an elevation in temperature of nanofluid. 
The increase in thermal conductivity raises temperature of nanofluid. 


AR WN 


The greatest increase in thermal conductivity is observed for bricks and least is observed for blades. 
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Blade shaped nanoparticles in blood show least velocity while brick shaped nanoparticles show highest 
velocity. 

The rise in the number of nanoparticles in the blood causes a decrease in velocity of nanofluid. 

The enhancement in the value of Grashof number brings a rise velocity of nanofluid. 

The enhancement in heat source parameter brings a rise in velocity of nanofluid. 


10. The enhancement in viscosity of nanofluid raises its velocity. 


11. The greatest increase in viscosity is observed for bricks and least is observed for blades. 


Thus, the shape of the nanoparticle can be selected based on the purpose of use and the respective 
thermodynamic properties can be controlled. The developed mathematical model has an avant-grade approach in 
the cure of cardiovascular diseases. 


6. 


Appendix 


The thermo physical properties of the blood are:- 


C, | 3594J/KgK 
p | 1063Kg/m* 

k | 0.492W/mK 
y | 0.18x10°1/K 


Table 1 
Shapes of nanoparticles; m | A 
Platelets 5./MBK1 
Blades 8.6 | 14.6 
Cylinders 4.9 | 13.5 
Bricks 3.7 | 1.9 
Table 2 


q4(z) = (legR(z) — logR,)*(4R(z) a Da(R(z)*(1 — 2R(2) + 2legR,) — RZ) 


+ VDa(A(z)*(—3 + 4legR(z) — 4logR.) + 4R(2)*R? — RF) 

+ A(z)a(A(z)*(1 — legR(z) + logR,) — 2R(z)*R? 

+ (1 + logR(z) — logR,)R7))/(16/(1 + Ab) (R(z)a(logR{z) — logR,) 
— VDa)) 


q2(z) = (logR(z) 


— logR,)" (2R¢ea (< (Re 4(—4+ 3lagR(z) — 3logR,) 

+ 4R(z)*(2 + logR(z) — logR,)R? 

- (4 + 7logR(z)+ 4(logR(z))” - (7 + SlogR(z) )logR, + 4(logR.)*)R?) 

A 2| ky, +(n—1)kp—(k,—ky)@ 
k,, + (m —1)k, + (n— 1)(k, - ke 

— R?) (R(2)*(6 — YogR(z) + 4logR(z))” + (9 - 8logR(z) )logR, 

+ 4(logR,)?) + 4R(z)?(—3 + (logR(z) — logR,)*)R? 


+ (é + 9lagR(z)+ 4logR(z)) - (9 5 BlogR(z) llogR, + s(looR.)") RS) 


(R(zi* 


+ VDa(12(R(z)*(9 — SlogR(z) + SlogR,) — 4R(z)*)(3+ 2logR(z) 
— 2logR,)(logR,)*(3+ 4logR(z) — 4logR,)(logR,)*) 

kyt (n- 1)ky _ (ky aks )e _\6 : 
+ ihe + @— 1k, + —1) (k, = | (—R(z) (27 — 46legR(z) 
+ 24(logR(z)) + 2(23 — 24logR(z)logR,+ 24iogR,)") 
+ 9R(z)*(7 — 4logR(z) + 4logR,)R2 — 9R(z)*(5+ 2logR(z)— 2logR,)R? 
+ (9+ BlogR(z) — 8logR,)R®)))) (768 (logR(z) — logR,))(R(2)a(logR(z) 
— logR,) — VDa)) 
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